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HIGH REYNOLDS NUMBER TESTS OF A 
C-141A AIRCRAFT SEMISPAN MODEL TO 
INVESTIGATE SHOCK-INDUCED SEPARATION 

By 

W. T. Blackerby and J. F. Cahill 


SUMMARY 


Tests of a semispan model of the USAF/Lockheed C-141 aircraft have been conducted 
in the Lockheed-Georgia high Reynolds number transonic wind tunnel (Compressible Flow 
Facility). Wing chordwise pressure distributions were measured over a matrix of Mach num- 
bers and angles of attack for which shock-induced separations are known to exist. The range 
of Reynolds numbers covered by these tests nearly spanned the gap between previously avail- 
able wind-tunnel and flight-test data. 

Results showed that use of the semispan test technique produced good correlation with 
the prior data at both ends of the Reynolds number range, but indicated strong sensitivity to 
details of the test setup. The severity of flow separation at low Reynolds numbers is shown 
to increase with spanwise distance from the centerline. The scale effect on that separation 
also increases with spanwise distance, resulting in a more uniform distribution of trai ling- 
edge pressure recovery at high Reynolds numbers. The rate of change of shock location and 
trai ling-edge pressure recovery with Reynolds number increases with increasing Reynolds 
number at outboard locations on the wing. The trends of these variations in severity of sep- 
aration are interpreted as an indication of spanwise contamination resulting from outboard 
flows in the boundary layer behind the shock. These data show that for the C-141 configuration no 
Reynolds number less than the full-scale value produces a satisfactory simulation of the flight test 
results. 


INTRODUCTION 


This reporf presenfs the results of tests of a semispan model of the USAF/Lockheed 
C-141 aircraft in the Lockheed-Georgia 20 X 28-inch high Reynolds number transonic wind 
tunnel (CFF). Wing pressure distributions were measured at three spanwise stations over a 
Reynolds number range from approximately 3 million to slightly over 20 million at Mach 
numbers from 0.75 to 0.85. 

The primary objective of this program is to develop data on the occurrence and the 
effects of shock-induced separation in the range of Reynolds numbers between those achiev- 
able in existing large transonic wind tunnels (AEDC 16-Foot Test TF-139, reference 2) 
and those experienced by the aircraft in flight (reference 1 ). Secondary objectives include 
a demonstration of the validity of semispan testing to obtain data of this type at high Rey- 
nolds numbers, an evaluation of whether a minimum Reynolds number can be identified as 
adequate for duplication of flight characteristics, and a correlation of data measured in this 
program with those from previous flight and wind tunnel testSo 

Several previous studies have been conducted on the shock-induced separation charac- 
teristics of the C-141 aircraft wing because a significant amount of data on sensitivity to 
scale effects was measured during the initial development of this aircrafto Results of those 
studies, presented in references 1 through 6, have confirmed the intimate relationship be- 
tween trailing-edge separation and shock location in these scale effects, the large influence 
of transition location, and a strong dependence on wing planform geometry (as distinguished 
from airfoil section effects). These facts are apparent in figure 1, which summarizes results 
from a large-scale panel model and from the complete aircraft configuration (reference 5). 
Tests of the panel model covered a range of Reynolds number from approximately 4 million 
to 28 million; but for the complete aircraft configuration, a gap was present in Reynolds 
number between the highest wind tunnel value of 8-1/2 million and the lowest flight value 
of 35 to 40 million. 

The Lockheed-Georgia CFF is designed to permit tests at Reynolds numbers up to ap- 
proximately 164 million per meter (55 million per foot)o A large semispan model of the 
C-141 aircraft designed for the CFF test section, 71.2 X 183 cm (20 X 28 in.), has a mean 
chord of 12.7 cm (5 in.) and, therefore, enables the acquisition of data at chord Reynolds 
numbers somewhat greater than 20 million. This program was undertaken to exploit the semi- 
span model technique and the capabilities of the CFF to expand current understanding of 
these transonic scaling phenomena. 
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SYMBOLS 
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Cp 
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MAC 
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wing span, cm (in.) 
pressure coefficient 

pressure coefficient at 30 percent chord 
lower-surface pressure coefficient 
trai ling-edge pressure coefficient 
upper-surface pressure coefficient 
freestream Mach number 
mean aerodynamic chord, cm (in.) 

local Mach number normal to the local element line of the wing 

2 2 

static pressure, N/m (Ib/ft ) 

ratio of static pressure to total pressure 

Reynolds number 

2 2 

wing area, m (ft ) 
nondimensional chordwise station 
shock location along the chord 
angle of attack, deg 

geometric angle of attack for the semispan model in the Lockheed-Georgia 
Compressible Flow Facility 

ratio of specific heats 

nondimensional spanwise station 

sweep angle of an element line on the wing 
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APPARATUS AND PROCEDURES 


Model 

The C-141 wing-fuselage semispan model used in this investigation was designed at 
the Lockheed-Georgia Company under contract NAS2-7081. This 0.0188-scale semispan 
model is shown installed in the 50.8 cm (20 in.) X 71 .2 cm (28 in.) working section of the 
Lockheed Compressible Flow Facility in figure 2. Model size is based on a good compromise 
between high Reynolds number requirements and consideration of model/wall interference 
and solid blockage. The model consists of a left wing, half fuselage, and gear fairing at- 
tached to a 0.51 cm (0.2 in.) thick plate as shown in the sketch of figure 3. The model 
plane of symmetry was displaced 3.8 cm (1.5 in.) from the tunnel floor to allow for the floor 
boundary layer, and the plate attached to the base provides reflection plane characteristics. 
The model span is 45.72 cm (18 in.), and the effective semispan/tunnel height ratio is 0.696 
with the 0.51 cm (1 .5 in.) displacement. The solid blockage is estimated to be 2.5 percent. 

Surface pressure orifices were installed on the wing at semispan stations of 0. 193, 
0.389 and 0.637 b/2. Orifices were provided on both upper and lower surfaces, and addi- 
tional orifices were located at intermediate stations to monitor spanwise effects. Orifice 
locations are detailed in table 1 . Spacing of the orifices was limited to 2.5 percent chord 
over most of the rear half of the airfoil at each of the three semispan stations. This ensures 
sufficient pressure data to accurately determine shock locations and local conditions in the 
vicinity of the shock. 


Instrumentation 

Measurements of the static pressures on the model surfaces were made using electron- 
ically actuated pressure scanning valves. The full-scale range of the quarter-percent accu- 
racy Statham transducers in the valves was selected to provide maximum accuracy for the 
wind tunnel conditions tested at model pressures of ± 34 dynes/cm^ (50 psi). CEC force 
balance pressure transducers were used in con|unction with CEC servo amplifiers to provide 
a precise measurement of the atmospheric pressure, stagnation pressure, and test section sta- 
tic pressure to 0.05 percent of the 172 dynes/cm^ (250 psi) capacity. These transducers 
allow determination of the test section Mach number to an accuracy of ± 0.002 at the high- 
est stagnation pressure. 

Angle of attack was measured with a calibrated potentiometer operated by the angle 
of attack drive mechanism. 

Raw pressure data were recorded on magnetic tape using the CFF high-speed data ac- 
quisition system. The data acquisition system consists of a Lockheed Electronics Company 
MAC- 16 computer and associated peripheral equipment. The raw data were reduced to co- 
efficient form with a CDC 1700 computer. 
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Test Facility 


The general arrangement of the Lockheed Compressible Flow Facility (CFF) is shown 

type control valve accurately maintains the settling chamber stagnation pressure at selected 
pressure less than or equal to the 172 dynes/cm^ (250 psia) maximum and at mass flow rates 
less than 1089 kg/sec (240 Ib/sec). The test section is 50.8 cm (20.0 in.) wide X 71.2 cm 
(28.0 in.) high X 183 cm (72.0 in.) long, and it is enclosed in a 3.7 m (12.0 ft) diameter 
plenum chamber. Two model support systems are available: a two-dimensional (side wall to 
side wall) system, and a three-dimensional wing or half-model system (bottom wall). The 
bottom-wall model mount system, used for semispan testing, has three variable-porosity walls 
and one fixed-porosity wall (bottom wall). Porosity is variable from zero to ten percent in 
the side and top walls by sliding two parallel plates with 0.635 cm (0.25 in.) diameter holes 
slanted 60 degrees from the vertical. The fixed-porosity floor has the same porosity pattern 
as that of the variable porosity walls, and the porosity can be reduced by manually sealing 
the holes. A more detailed description of facility is presented in reference 7. 


in figure 4. The tunnel is of the blow-down t^ye, exhausting directly to the atmosphere 
The air storage capability is 368 m^ (13,000 ft^) at 413 dynes/cm^ (600 psia) . A sleeve 
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TESTS 


General Test G>nd!t!ons 

Unit Reynolds number for the test ranged from a low of about 19.7 million per meter 
(6 million per foot) to o maximum of about 190 million per meter (58 million per foot). 

This corresponds to a variation in model Reynolds numbers from about 2.5 million to 24 mil- 
lion/MAC. The Reynolds numbers were selected to provide wing shock/boundary layer char- 
acteristics which close the gap between model and full-scale results. Reynolds numbers for 
previous wind tunnel tests on the complete model were limited to 8.5 mil lion/MACand flight 
test Reynolds numbers range from 35 to 70 mil lion/MAC. Mach numbers for the current 
program were 0.75, 0.8, 0.825, and 0.85; these were selected, together with model geo- 
metric angles of attack of 1^, 2®, 3®, and 4°, to correspond to those flight conditions 
where scale effects on wing shock characteristics are prominent. A summary of the current 
test program Is given in table 2. 

The model was tested smooth and with transition strips of Ballotini glass beads. The 
Ballotini strips were placed on upper and lower wing surfaces and on the fuselage nose. A 
particle height of 0.0069 cm (0.0027 in.) was determined for the low-Reynolds-number test 
condition, according to the criteria of reference 8, and was used throughout the fixed tran- 
sition portion of the test. The roughness strips were placed 0.76 cm (0.3 in.) aft of the 
wing leading edge and 1 .9 cm (0.75 in.) aft of the fuselage nose to simulate the location 
used in reference 2. The Ballotini strips were 0.13 cm (0.05 in.) wide and were set with 
Eastman 910 adhesive. 


Angle-of-Attack Correlation 

Difficulties in measuring angle of attack accurately during maneuvers in flight tests, 
differences in aeroelastic twist of the various models considered and of the flight test aircraft, 
differences in model configurations, and the lack of a precise method for calculating wall- 
induced flow angles in perforated wall transonic wind tunnels combine to produce a require- 
ment for definition of an effective angle of attack which will assure that data from all avail- 
able sources are being compared for equivalent flow conditions. The basis for the concept used 
in this study is illustrated by the data from another test program shown in figure 5. "^his figure 
shows a comparison of data at one test condition from a wind tunngl test at 7.4X 10^ Reynolds 
number with data from flight tests at a Reynolds number of 82 X 10 . Because of the difference in 
Reynolds number, a significant difference in the degree of trailing-edge separation occurs, and 
the shock location changes by nine percent of the wing chord. In spite of these differences in 
shock location and in local velocities over the after portion of the wing, the remainder of the 
pressure distribution is in nearly exact agreement. On the other hand, variations in angle of 
attack (as shown by plots in the Appendix) cause large differences in the entire level of the 
velocity distribution. The concept used in defining the effective angle of attack, therefore, 
depends on the following assumptions, which are well supported by data from this and several 
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prior studies: 

1 «, At a given Mach number, only angle of attack causes a change in pressure dis- 
tribution forv/ard of the shocko With fully developed transonic flow, changes in 
circulation due to trai ling-edge separation cause no significant changes in pres- 
sure distribution except in the subsonic flow behind the shock. 

2o Since the shape of the pressure distribution remains constant at any fixed Mach 
number and angle of attack, specification of pressure coefficient at any single 
point on the v/ing upper surface defines an angle of attack® 

The method of implementing this concept is illustrated by figure 6. In the upper part 
of this figure, the value of the pressure coefficient at 30 percent chord, ^P^ 3 / plotted 
against angle of attack for a Mach number of 0.825. In the lower portion of the figure, 
shock locations from previous wind tunnel and flight tests are plotted against their accom- 
panying values of Cp Values of shock location selected from these plots at values of 
Cp^ 2 corresponding to the angle of attack for the CFF tests then enables a proper compari- 
son'of data from these other sources with data from the CFF. This determination of equiva- 
lent angle of attack was accomplished for each spanwise station individually and, therefore, 
compensated for differences in aeroelastic twist among the several data sources. 

Model and Tunnel Configuration 

Prior to the start of the basic test program, a series of preparatory test runs were made 
to verify the proper model and tunnel configuration. The rationale for selection of the final 
test configuration v/as based on subfective Judgments made after comparing current semispan 
test data with previous complete model results. A good match of pressure data from AEDC 
Test TF-139 (reference 2 ) on the wing upper surface back to the shock was the basic criterion 
used to define an equivalent geometric angle of attack® Consideration was then given to 
matching shock location at the equivalent conditions. Since correlation of the lower surface 
pressures is somev/hat limited by configuration differences, no significance can be attached to 
either local lower surface pressures or lift coefficients for this analysis® A brief discussion of 
these differences is given in the next section. 

The initial tunnel/model configuration was as follows: 

o Model and model support fang as shown in figure 3® The model support between 
the boundary layer diverter and the tunnel floor having a simple rectangular 
cross-section. 

o CFF floor porosity of six percent, other walls having a porosity of four percent. 

The floor porosity was selected as a starting point on the assumption that it represented 
a reasonable compromise between a 100 percent reflection plane with a thick boundary lay- 
er and a partial (but adequate) reflection plane from which a substantial portion of the bound- 
ary layer bleeds through the partial porosity. A four-percent porosity on the other three walls 
v/as selected based on experience gained from recent testing of two-dimensional airfoils in the 
CFF (reference 9). The rectangular cross-section of the model support tang was assumed to 
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be satisfactory because of shielding provided by the boundary layer diverter. Results ob- 
tained for this initial configuration indicated that the partial floor porosity was undesirable, 
and a quick check was obtained by testing with tape placed over the floor holes. The floor 
was then sealed firmly by filling all the holes with epoxy, and the balance of fhe testing 
was accomplished with the floor sealed in this manner. 

Figure 7 shows preliminary wall porosity effects obtained after sealing the floor. 
Semispan data at two degrees geometric angle are compared with complete model data from 
AEDC TF-139 at zero degree geometric angles This two degree difference was selected to 
approximately match the flow condition on the two models and is probably due to the difference 
in wall-induced angularity between the two tests discussed previously. These data indicate 
expected increases in local velocities over the airfoils accompanied by rearward shock move- 
ment as the wall porosity is gradually decreased. At the extremely low porosity of one per- 
cent, noticeable velocity increases have occurred, and the trail ing-edge pressure recovery 
has deteriorated substantially, especially on the two outboard stations. The aft shock move- 
ment with lower porosity improves the agreement of semispan and full-span pressure distribu- 
tions and implies use of a rather low porosity for the test. At this point in the preliminary 
testing, a value of two percent was chosen while further checks on the overall tunnel/model 
configuration were being made. 

Figure 8 summarizes the effects on the chordwise pressure distributions for the mafor 
steps taken to finalize the configuration. Comparisons of CFF results at two percent poro- 
sity are made against data from AEDC TF-139, again at equivalent angle of attack. Agree- 
ment of data for the initial configuration with full-span results is confined to the most for- 
ward portion for the two inboard stations with some improvement outboard. Shock locations 
for the initial configuration are six to ten percent forward of the full-span case. The sub- 
stantial improvement in the correlation gained by sealing the floor can be clearly seen in 
these comparisons. A smaller improvement was obtained by adding a fairing to the model 
support tang below the diverter plate. This fairing changed the cross-section of the tang 
from rectangular to that of an airfoil shape and significantly increased its fineness ratio. 
Because of the improvements due to these two changes, the test configuration was updated 
to include both the sealed floor and the support tang fairing. 

Before concluding the preparatory testing, a final check on porosity was made with 
this test configuration using values of four, three, and two percent for wall porosity. The 
results are shown in Figure 9. in this case, a porosity of two percent produced the substan- 
tia! loss in trailing-edge pressure recovery seen previously at one percent. The reason for 
this cannot be substantiated but is associated with large reductions in the drag of the overall 
model/support configuration due to the addition of the tang fairingso Based on these data, 
selection of three percent as the correct porosity was made, and the main program was run 
using this value. The tang fairing used for these test runs and all subsequent testing was 
slightly thicker than that for the data discussed in the previous paragraph and produced 
slightly different test results. 
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Correlation of Pressure Distribution with Full-Span Data 

As mentioned previously, lower surface pressures measured on the semispan model do 
not correlate with data measured on the complete model in AEDC TF-139c This difference, 
together with other observations regarding correlation of the CFF pressure data, will be 
discussed here* 

Figure 10 shows an overall comparison of pressure results from three sources: current 
CFF semispan test, AEDC TF-139 test, and flight test. These data have approximately the 
same pressure coefficient at 30 percent chord on the Inboard part of the wing^ This assures 
an equivalent angle of attack inboard. Any spanwise variation in local angle of attack will 
result in further discrepancies outboard. Figure 10 reveals that this does, in fact, happen. 
The data for the outboard station indicates that this section on the aircraft is obviously at an 
effective lower angle of attack than the same section on the AEDC and CFF models. This 
is due purely to differences in aeroelastic twist on the wind tunnel model wing compared 
with the full-scale aircraft wing. Typical values of aeroelastic twist calculated for condi- 
tions encountered in the flight test program are 1^1/2 to 2 degrees at the wing tip. 

Two significant differences in the wind tunnel models can be identified which at least 
partially explain the kind of flow differences seen on the inboard lower surfaces in figure 10: 

o The model used in AEDC TF-139 consisted of a C-141 wing mounted on a C-5 
model fuselage. 

o The semispan model did not include pylons and nacelles. 

The exchange of a C-5 fuselage and wheel pod fairing for those of a C-141 could 
produce noticeable flow changes on the inboard lower surface of the wing. The differences 
in pressure distribution on the lower surface of the inboard station, aft of 30 percent chord, 
between the flight data and AEDC TF-139 test data indicate the effect of the fuselage/ 
wheel pod exchange. The peak pressure for the flight case is shifted further aft with respect 
to the AEDC data, and the aft loading is increased. The semispan data, where the C-141 
fuselage/wheel pod is also used, tends to follow the pattern of the flight data. The large 
discrepancy forward on the section is obviously not being caused by the fuselage/wheel pod 
exchange. Some of this difference can be attributed to the absence of pylons and nacelles 
on the semispan model. Previously measured data on a similar configuration indicate that 
approximately half of the increment between the complete and semispan model results could 
be due to pylon and nacelles. 


RESULTS 

The basic test results are shown in Appendix A in the form of computer-generated 
pressure distribution plots. Generally, each figure presents pressure distributions for a series 
of angles of attack at fixed values of Mach number, Reynolds number, and spanwise location. 
Typical developments of flow separation as the angle of attack or Mach number increase can - 
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be observed In a gross sense simply by comparing these pressure distribution plots. For ex- 
ample, Appendix figures 52A to 54A, showing data measured at a Reynolds number of 15 
X 10^ and a Mach number of 0«825, indicate a progressive deterioration of trailing-edge 
pressure recovery as the angle of attack is increased, with an accompanying forward move- 
ment of shock location. These effects Increase in severity with increasing spanwise distance 
from the centerline. Figures 25A to 27A, for a Mach number of 0.80, show substantially 
smaller separation effects; and figures 76A to 78A, fora Mach number of 0.85, indicate 
somewhat greater separation. Details of the scale effects displayed by these data will be 
presented in the body of the report. 

Definition of Shock Location 


Figure 11 illustrates the definition of shock location used in this report. The abrupt 
pressure rise through the shock is fitted with a straight line which is extrapolated, if neces- 
sary, to a pressure coefficient corresponding to a local Mach number of 1.0 considered nor- 
mal to local element lines of the wing, 
ing equation: 

2 


r/1 


Cp = 


YM^ 


This pressure coefficient is defined from the follow- 

Y 

- 1 


+ m2 2 

+ — cos 


1 


+ — M|_n 


where is the sweep angle of each element line on the wing (constant percent chord), 
Mln '5 the Mach number corresponding to the local velocity on the wing surface normal to 
the element line (set equal to 1.0 for this purpose). (Other terms have their conventional 
meanings and are defined in the Symbols Section.) The chordwise location of the intersec- 
tion of the extrapolated shock line with the local Mach number 1 .0 line is defined as the 
shock location, X/C^i^. This concept of considering Mach numbers normal to local element 
lines is used at this point only as a convention for defining shock location and is completely 
valid for this purpose when applied in a consistent manner. If Mach numbers normal to iso- 
bars were used instead, a foundation would exist also for a quantitative prediction of shock 
characteristics. In the absence of sufficient data to define isobar sweep, the sweep of ele- 
ment lines is frequently used as a reasonable approximation for those quantitative purposes. 

Basic Scale Effects 


Using the angle-of-attack correlation concept described previously, data from pre- 
vious investigations have been combined with those from the present tests to show scale ef- 
fects on shock location and trailing-edge pressure coefficient over the full range of avail- 
able data. These data are presented in figures 12, 13, and 14 for each of the Mach num- 
ber/angle -of-at tack combinations for which strong shocks were observed. Because of the 
unsteady nature of the separated flows which are being investigated, a modest amount of 
scatter appears in all of the data shown. However, this scatter is not so large as to obscure 
the trends of the data, and good correlation generally exists among the data from the three 
sources available. In several instances (spanwise station .637, figure 12(a) is a good ex- 
ample), two distinct variations in shock location are shown for low Reynolds numbers with 
the transition fixed. It is believed that the differences shown can be attributed to the fact 
that, in some of the test runs, the boundary layer transition occurred at the grit strip, while 
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in ofhers, fhe transifion pornf moved farfher forward. The curves drawn on these plots are 
faired through the points showing the more forward shock locations, since the higher Rey- 
nolds number results (where the nature I -transition and the fixed-transition data agree) 
should experience the more forward transition location. A similar variation is shown by 
drag data presented in reference 9. A systematic disparity between the shock locations ob- 
served in flight and those which would be expected from the semispan results is shown for 
the two inboard stations at low angles of attack for a Mach number of 0o8. The source of 
this apparent discrepancy has not been identified* At a Mach number of 0.85, where 
flight data were measured at Reynolds numbers of 45 and 60 million, the trailing-edge 
pressure recovery is poorer, and the shock location is farther forward at the higher of these 
Reynolds numbers. Since this trend is in opposition to logical expectations, the curves are 
faired between the two values. 

Low Reynolds number data from tests with natural transition indicate a tendency, 
shown in several previous studies (references 3 and 10),to provide a better simulation of the 
flight results than do those with the transition artificially fixed. This tendency is demon- 
strated also by the data shown in figures 12 to 14 (except figure 12a where the flight shock 
locations appear suspiciously far forward). The difference between the fixed transition and 
the natural transition data is significantly less pronounced for the most inboard wing station, 
where the flow can be contaminated by boundary layer flow from the fuselage. In the ab- 
sence of sufficiently high Reynolds number test facilities, testing at low Reynolds numbers 
with natural transition still appears to provide a useful means to improve the accuracy of 
full-scale simulation. Unfortunately, quantitative methods to provide a precise simula- 
tion of a specific high Reynolds number condition have not yet been demonstrated. 

As discussed previously, the model used in the present tests did not include pylons 
and nacelles, while all data available from previous testing were obtained with the pylons 
and nacelles installed- The good correlation shown at both high and low Reynolds number 
in figures 12 to 14 indicates that the absence of the pylons and nacelles has little, if any, 
influence on upper surface flows in spite of the obvious effects on lower surface pressure 
distribution. An initial criterion in C-141 development specified that the pylon leading 
edge should intersect the wing surface aft of the attachment line. Therefore, it is logical 
that the upper surface boundary layer is not disturbed by the pylons- 

Evidence of Spanwise Contamination 

It was previously pointed out that the tro*! ling-iedge separation increases as the mea- 
surement station moves outboard. This fact is strongly reflected by the much more forward 
shock location for outboard stations at low Reynolds numbers in figures 12 through 14. 

These figures also show a more pronounced scale effect on separation and shock location 
for the outboard stations. Therefore, the spanwise variation in shock location is very much 
smaller at high Reynolds numbers* 

In a number of instances, especially at the outboard stations, the data in figures 12 
through 14 show an increasing slope in the scale effect on shock location at Reynolds num- 
bers above 15 to 20 million. Figures 15 and 16 combine data from several of those figures to 
show how this slope increase develops with increasing Mach number or angle of attack. 
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Figure 15 shows data for several angles of attack at a Mach number of 0.825. At 
one degree angle of attack, the scale effect curves have a nearly constant slope for station 
.637, culminating in essentially unseparated flow at the highest Reynolds number. For the 
two more inboard stations, the pressure recovery reaches typical unseparated values at low- 
er Reynolds numbers, and the shock locations show a plateau for Reynolds numbers of ap- 
proximately 30 to 60 million. 

At two degrees angle of attack, the pressure recovery still shows attached flow at 
station .637 and 40 million Reynolds number; but the separation increases rapidly as the 
Reynolds number is decreased, causing the shock location to move rapidly forward for Rey- 
nolds numbers down to approximately 10 million. The shock location curve shows evidence 
of approaching a plateau at station .389, and a distinct plateau still exists at station .193. 

At higher angles of attack, no flight data are available for station .637, but the 
data for station .389 indicate increasing slopes of the aft movement of shock location for 
Reynolds numbers above approximately 15 million. A plateau is still apparent at station 
.193 for high Reynolds numbers, even for the highest angle of attack. 

Figure 16 shows scale effect data at an angle of attack of two degrees for several 
Mach numbers. At station .637, essentially attached flow is indicated for the highest 
Reynolds number and for all Mach numbers. The separation again increases rapidly at .825 
and .85 Mach numbers as the Reynolds number decreases, resulting in an abrupt variation 
in shock location with Reynolds number. At station .389, the separation for the highest 
Reynolds numbers is more severe than at station .637, and no abrupt movement in shock 
location occurs. At station .193, the separation is never very severe, and a normal aft 
movement of shock location with increase in Mach number is shown for the highest Reynolds 
number. 

These variations would seem to indicate the presence of spanwise contamination ef- 
fects, causing separations to develop more rapidly than anticipated at outboard stations be- 
cause of outboard flows along the trailing edge. At inboard stations, the separation is pri- 
marily a function of the local section and planform properties. At outboard stations, those 
normal separation tendencies are amplified by additional influences contributed by the 
spanwise flow. For the highest Reynolds numbers, of course, the spanwise contamination 
effects are minimal because no strong separations are present on the inboard portions of the 
wing. 


Shock locations and trailing-edge pressure coefficients for two degree angle of attack 
at .825 Mach number are shown cross-plotted against spanwise station in figure 17(a). At 
a Reynolds number of 10 million, the trailing-edge pressure recovery deteriorates rapidly 
from inboard to outboard stations, and the shock location moves forward approximately 20 
percent of the chord. At 40 million Reynolds number, the trailing-edge pressure coeffi- 
cient is near 0.2 at all spanwise stations, and the shock location changes by only three 
percent chord. The pressure recovery is slightly worse at station .389 than at either the 
more inboard or the more outboard station in this case, possibly as a result of the break in 
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the trailing-edge sweep angle which is quite near this station. 

In figure 17(b), similar data are shown for a Mach number of 0.85. Again, a rapid 
increase in separation and forward shock movement occur for increasing outboard location 
at 10 million Reynolds number. At higher Reynolds numbers, the separation is obviously 
most severe at station .389, with significant improvement at the most outboard station. 

This fact does not contradict previous statements regarding spanwise contamination, since 
the effects of increased Reynolds number in improving the resistance to local separation 
tendencies can overpower the effects of spanwise flow. 

Figure 18 shows the pattern of surface flows on a large-chord panel model having no 
taper and no spanwise variation in airfoil section or twist. The airfoil section of this model 
was the same as that at wing station .389 on the C-141A wing, and three-dimensional 
streamline-shaped plates were attached to the root and tip of the model to eliminate end 
effects. Near the root, the flow is smooth and continuous through the shock (indicated by 
the abrupt change in flow direction at approximately 0.75 chord). A strong spanwise flow 
appears in the boundary layer behind the shock. At more outboard spanwise stations, the 
shock moves forward on the chord, and a distinct discontinuity appears in the boundary lay- 
er flow at the shock. Although not conclusive, this would appear to offer additional evi- 
dence of the influence of spanwise flows on the severity of trailing-edge separations and, 
therefore, on shock location. 


Correlation of Shock Location and Trailing Edge Separation 

Previous references have been made to the direct correlation of trailing-edge separa- 
tion with shock location, and this is illustrated by the similarity in the shape of sea I e-effect 
curves for these two parameters. In figure 19, shock location is plotted against trailing- 
edge pressure coefficient using data from the various sources available and showing that, in 
a quantitative sense, this correlation is quite rigorous. These data lend additional weight 
to the argument that the kind of scale effects considered here are predominantly a function 
of trailing-edge separation effects rather than discrete shock -boundary layer interaction 
phenomena in the immediate vicinity of the shock. In other instances (with less severe aft 
pressure gradients, for example), the effects at the shock might be dominant, but no data 
from the present study indicate that kind of phenomenon. 

The data shown in figure 20 show further substantiation of the strong dependence of 
shock location on trailing-edge separation phenomena and illustrate the mechanism by which 
the shock is caused to move. In figure 20(a), data are shown from reference 2 with, and 
without, vortex generators at 70 percent chord on a C-141 wing model. With no vortex 
generators, a small region of separated flow is apparent aft of 85 or 90 percent chord. 
Adding vortex generators aft of the shock eliminated that separation and caused the shock to 
move rearward. The direct cause of the shock movement is the change in the pressure distri- 
bution from the foot of the shock aft resulting from elimination of the trailing-edge separa- 
tion. The shock acts to produce the change in local Mach number from supersonic to sub- 
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sonic. Since the pressure rise through the shock is a function of the supersonic Mach num- 
ber at which it occurs, it must position itself in such a way that its pressure rise is compat- 
ible with both the upstream local supersonic Mach number and the downstream pressure dis- 
tribution. Any change which causes the downstream pressures to be more negative, as 
shown here due to eliminating the trailing-edge separation, requires the shock to occur at 
a higher local Mach number and, therefore, to move aft. 

In figure 20(b), a similar change In shock location is shown to be caused by changes 
in downstream pressures resulting from aileron deflection. Upward (negative) deflections of 
the aileron cause the pressure on the rear of the wing to be more positive and, therefore, 
require the shock to move forward. The fact that the downstream pressures comprise the 
driving factor and that shock location is the result (rather than the inverse) is demonstrated 
by the fact that the downstream pressures are predicted quite well by a Karman-Tsien trans- 
formation of pressures measured for the same aileron deflections at low Mach numbers for 
which no shock was present. 

In the case of the scale effects under consideration here, no change in airfoil geome- 
try is involved; as the pressure changes downstream of the shock are uniquely controlled by 
trailing-edge separation, it is possible to correlate shock location against trailing-edge 
pressure recovery. 


Comparison with Panel Model Data 

Pressure distribution data measured at station .389 at a Reynolds number of 20 million 
are compared in figure 21 with data measured at three stations near the mid-semispan of the 
panel model of reference 5. Excellent agreement of the pressure distributions is shown over 
most of the chord. Immediately forward of the shock, a rather rapid increase in velocity is 
shown in the data from the present investigation. The velocity variation in this small region 
is accurately simulated by the panel data at station 40. At the other two panel model sta- 
tions, this small area of increased velocity is less noticeable. 

Figure 22 shows scale effects on shock location from the present tests and from the 
panel model test for two angles of attack at M = .825. At the highest Reynolds number, or 
at. the lowest Reynolds number with natural transition, the shock locations from these two 
sources agree very well. However, at lower Reynolds numbers with the transition fixed, the 
separation is more severe in the complete configuration data, and the scale effect curves 
are, therefore, significantly steeper. 

This failure of the panel model data to simulate scale effects on the complete config- 
uration, when viewed in conjunction with other previously discussed effects, would appear 
to result from its inadequate simulation of the complete configuration geometry (either in- 
adequate span for matching span wise contamination or the absence of the local disturbance 
produced by the break in trailing-edge sweep angle). The surface flow pattern shown in 
figure 18 indicates that, if measurements were available at more outboard stations on the 
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panel, they would show more pronounced separations and certainly more forward shock 
locations. Caution is obviously required in attempting to predict complete wing scale 
effects on the basis of data from simplified models such as this panel or airfoil sections. 



CONCLUSIONS 


1. Good correlation of data from several sources indicates that these semispan tests produced 
a good simulation of flow conditions existing on the complete configuration. However, pre- 
paratory testing to establish the final model/tunnel configuration showed a strong sensitivity 

to configuration details and indicated that a large boundary layer diverter plate cannot be relied 
upon to shield the model from inboard disturbances. 

2. Results from this program provide a good bridge between previously available low-scale wind 
tunnel data and full-scale flight data. Abrupt increases in the rate of variation of shock 
location with Reynolds number are shown for high Reynolds numbers, especially on outboard 
portions of the wing span. 

3. At low Reynolds numbers, the severity of trailing-edge separation increases with increase in 
spanwise distance from the centerline, apparently as a result of spanwise flows in the region 
behind the shock. At higher Reynolds numbers, this spanwise contamination is significantly 
reduced. 

4. For the Cl 41 configuration, no Reynolds number less than the full-scale flight value 
provides a satisfactory simulation of the full-scale flows. Furthermore, the abrupt increase 
in the rate of variation of shock location with Reynolds number increases the difficulty of 
reliable extrapolation of low-scale measurements. 

5. Because of the spanwise contamination effects suggested by the data from these tests, data 
from simplified models such as constant-chord panels or two-dimensional models cannot be 
relied upon to produce an accurate indication of scale effects involving shock-induced 
separation. 
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TABLE 1 


ORIFICE LOCATIONS 


STATION 

.193 

.298 

.343 

.389 

.422 

.488 

.637 

X/C 

UPPER 

LOWER 

UPPER 

UPPER 

UPPER 

LOWER 

UPPER 

UPPER 

UPPER 

LOWER 

0 

o 

o 



o 

o 



o 

0 

.025 

o 

o 



o 

o 



o 

o 

.05 

o 

o 



o 

o 



o 

o 

.075 

o 

0 



o 

o 



o 

o 

.1 

o 

o 



o 

o 



o 

o 

.125 

o 




o 




o 


.15 

o 

o 



o 

o 



o 

o 

.175 

o 




o 




o 


.2 

o 

o 



o 

o 



o 

o 

.25 

o 




o 




o 


.3 

o 

o 



o 

o 



o 

o 

.35 

o 




o 




o 


.4 

o 

o 



o 

o 



o 

o 

.45 

o 




o 




o 


.475 

o 




o 




o 


.5 

o 

o 



o 




o 

0 

.525 

o 




o 




o 


.55 

o 




o 




o 


.575 

o 




o 




o 


.6 

o 

o 



o 

o 



o 

o 

.625 

o 




o 




o 


.65 





o 




o 


.675 

o 




o 




o 


.7 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

.725 

o 




o 




o 


.75 

o 




o 




o 


.775 

o 




o 




o 


.8 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

.85 

o 




o 




o 


.9 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

.95 

o 




o 




o 


1.0 

o 

o 



o 

o 



o 

® 1 
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TABLE 2. PROGRAM SUMMARY 


WING 

ReXlO^ 

MACH 

RUN NUMBER 

@ INDICA 

TED AI^LE 

OF ATTAC 

GRIT 

MAC 

NO. 

1° 

2© 

? 

4° 

FREE 

2.7 

.75 

394, 392* 393 

396, 395* 

397 

398 



.8 

386 

387 

389, 388* 

391 , 390* 



.825 

412(409-411)* 

414,413* 

415 

416 



.85 

380 

381 

427, 426* 

385, 384* 


4.4 

.825 

374 

375 

376 

378, 377* 



.85 

417 

418 

419 

420 


6.0 

.75 

201 

202 

203 

224, 222* 







223* 



.8 

205 

206 

207 

208 



,825 

211,209,210* 

212, 402 






401, 399* 400* 

425 

213, 403 

214, 404 


• 

.85 

216,215* 

217 

218 

221,219* 




405 

406 

407 

220* 408 


8.2 

.75 

226, 225* 

227 

228 

229 



.8 

230 

231 

232 

233 



.825 

234, 421 

235, 422 

236, 423 

237, 424 



.85 

240, 238* 239* 

242 

243 

244 


15 

.75 

246, 245* 

247 

248 

249 



.8 

250 

251 

252 

253 



.825 

254 

255 

256 

258 



.85 

259 

261 

262 

263 


20 

.75 

265 , 264* 

266 

267 

268 



.8 

273 

274 

275 

276 



.825 

269 

270 

271 

272 



.85 

277 

278 

279 

280 


24 

.825 


281 



FIXED 

2.7 

.8 

444 

445 

446 

447 



.825 

432(429-431)* 

433 , 323 

434 

436, 435* 



.85 

440(437-439)* 

441 

442 

443 


6 

.75 

334 

336, 335* 

337 

338 



.8 

340, 339* 

341 

342 

343 




449 

450 

451 

452 



.825 

330 

331,325 

332 

333 





324* 





.85 

344 

345 

346 

347 




453 

455, 454* 

456 

457 


8.2 

.75 

352 

353 

354 

355 



.8 

357, 356* 

358 

359 

360 



.825 

348 

349, 327 

350 

351 





326* 





.85 

362 

363 

364 

365 


15 

.825 

366 

367, 328 

368 

369 



.85 

370 

371 

372 

373 


20 

.825 


428, 329 


458 


REMARKS 


(a=0” RUN 200) 
(a=0°, RUN 204) 


GRIT ERODED ON 
RUNS 339-343 


GRIT ERODED ON 
RUNS 344-347 


NOTE: ASTERISKS DENOTE ADDITIONAL MACH NUMBER RUNS. 















EL FLOOR 


10. 16 cm 
(4 in.) 


7,62 cm 
(3 In.) 




7.62 cm 
(3 in.) 


C-141A Semispan Model 











I DODO 


DATA 

REYNOLDS NO. 

MACH NO. 

ALPHA 

POROSITY 

CFF 

7.5 X IOVmAC 

0.8 

20 

6% 

CFF 

7.5 X IoVmAC 

0.8 

2° 

4% 

CFF 

7.5 X IOVmAC 

0.8 

2° 

2% 

CFF 

7.5 X lOVMAC 

0.8 

' 2® 

1% 

TF-139 

8.5 X IOVmAC 

0.8 

0° 







DATA REYNOLDS NO. MACH NO. ALPHA 


o 

CFF 

7.5 X IoVmAC 

0.8 

2° 

INITIAL DATA 

□ 

CFF 

7.5 X lOVMAC 

0.8 

20 

FLOOR SEALED , 

A 

CFF 

7.5 X lOVMAC 

0.8 

20 

FLOOR SEALED & 
SUPPORT TANG FAIRED 


TF-139 

8.5 X IOVmAC 

0.8 

0° 




0 0.2 0.4 0.6 0.8 1.0 

X/C 


Figure 8. Model/Tunnel ConfiguraHon Effects 








DATA REYNOLDS NO. MACH NO. ALPHA 


O CFF 7.5 X IoVmAC 

□ CFF 7.5X10VMAC 

A CFF 7.5 X 106 /MAC 

— . TF-139 8.5 X 106 /MAC 


0.8 

0.8 


2° INITIAL DATA 

2® FLOOR SEALED 

79 FLOOR SEALED & 

SUPPORT TANG FAIRED 

00 


0.8 





/ 

/ 

/ 

/ 


DATA 

REYNOLDS NO. 

MACH NO. 

ALPHA 

POROSITY / 

O CFF 

7.5 X IOVmAC 

0.8 

2° 

4% / 

□ CFF 

7.5 X IOVMAC 

0.8 

20 

3% 

A CFF 

7.5 X IOVmAC 

0.8 

2© 

2% 

TF-139 

8.5 X IOVmAC 

0.8 

0° 




x/c 


Figure 9. Porosity Effects, Floor Sealed and Support Tang Faired 
(a) n = .193 
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DATA 

REYNOLDS NO. 

MACH NO. 

ALPHA 

o 

CFF 

7.5 X IoVmAC 

0.8 

20 

□ 

CFF 

7.5 X K)6/MAC 

0.8 

20 

A 

CFF 

7.5X10VmAC 

0.8 

20 

— 

TF-139 

8.5 X lOVMAC 

0.8 

0° 


POROSITY 

4% 

3% 

2% 










.8 i:o 


0 .2 


4 .6 

X/C 


.8 1.0 


L -.6 

X/C 

e Data with Full Span Data, 
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Figure 16. 


Composite Summary of Flight and Wind Tunnel Scale Effects at 2 


o 


Angle of Attack 
















Oi 


M = .8 

a = 0® 



BEHIND SHOCK, REFERENCE 2 


M = 0.825 



(b) EFFECT OF AILERON DEFLECTION ON 
SHOCK LOCATION, REFERENCE 11 


Figure 20. Downstream Pressure Controls Shock Location 





Figure 21. Comparison wi|-h Panel Model, M 







□ SMOOTH 
O TRANS FIXED 
^ AEDC TF-139, REF. 2 
A FLIGHT, REF. 1 


O lU 

Rg - IoVmAC 

4U 

80 X 10' 

Figure 22. Concluded 







APPENDIX A 

MACHINE-PLOTTED PRESSURE DISTRIBUTION DATA 



The plotting sequence follows the program summary of table 2 in the body of this re- 
port. In general, each figure contains plots for the angle-of-attack variation at a span- 
wise station for each measured condition. Data for all three spanwise stations are plotted 
in sequence, and the figures are numbered consecutively. A summary of figure numbers is 
given in table lA. 



TABLE lA. FIGURE INDEX FOR PLOTTED DATA 



FREE TRANSITION 


Re X 10“6| SPAN STATION 


ON MAC 



.389 .637 


Fig. Fig. 
2A 3A 
5A 6A 

I 

8A 9A 
llA 12A 
14A 15A 


16A 17A 18A 

19A 20A 21 A 
22A 23A 24A 
25A 26A 27A 
28A 29A 30A 
31 A 32A 33A 



35A 36A 
38A 39A 
41 A 42A 


43A 44A 45A 
46A 47A 48A 
49A 50A 51 A 
52A 53A 54A 
55A 56A 57A 
58A 59A 60A 
61 A 62A 63A 
■ 65A 66A 


68A 69A 
71 A 72 A 
74A 75A 
77A 78A 
80A 81A 


MACH 

NO. 


.75 1®-4® 

.75 1°-4® 

.8 1°-4® 

.8 1®-4° 

.8 1®-4® 

.8 1°-4° 

.825 1°-4° 
.825 1°-4° 
.825 1°-4® 
.825 1®-4° 


FIXED TRANSITION 


R^ X 10-6 SPAN STATION 
ON MAC 1.193 ( .389 

Fig. I Fig. 
83A 
85A 86A 
88A 89A 
91 A 92A 
94A 95A 
97A 98A 
100A 101 A 
103A 104A 
106A 107A 
109A 110A 
112A 113A 
115A 116A 
118A 119A 
121A 122A 
124A 125A 
127A 128A 
15.0 130A 131A 
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